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Shedding light on astrophysical and particle physics 
origins of cosmic-rays can lead to breakthroughs in our 
understanding of the fundamental laws ruling the uni- 
verse. The goal of this paper is to introduce a new model 
independent approach for combining observations from 
different cosmic-rays experiments. We show that our re- 
sults lead to predictions relevant for future observations 
and to constraints useful to guide model builders. The 
data recently collected by PAMELA pQ indicate that 
there is a positron excess in the cosmic ray (CR) en- 
ergy spectrum above 10 GeV. The rising behavior ob- 
served by PAMELA does not fit previous estimates of 
the CR formation and propagation implying the possible 
existence of a direct excess of CR positrons of unknown 
origins. Interestingly PAMELA'S data show a clear fea- 
ture of such a positron excess but no excess in the anti- 
protons. While ATIC [2j and PPB-BETS [3j reported 
unexpected structure in the all-electron spectrum in the 
range 100 GeV- 1 TeV, the picture has changed with the 
higher-statistics measurements by FERMI-LAT [3] and 
HESS [5], leading to a possible slight additional unknown 
component in the CR flux over and above the specific 
Moskalenko and Strong model prediction [6j [7] which as- 
sumes a single-power-law injection spectrum of e ± cosmic 
rays. These interesting features have drawn much atten- 
tion, and many explanations have been proposed: For 
example, these excesses could be due to an inadequate ac- 
count of the CR background in previous modeling; The 
presence of new astrophysical sources; They could also 
originate from annihilations and/or decays of dark mat- 
ter. We refer to [8] for a recent review. Whatever the 
origin of these excesses might be, we show that we can 
derive novel constraints able to shed light on the physical 
nature of their source and/or propagation. We start by 
writing the observed flux of electrons and positrons as 
the sum of two contributions: A background component, 
<j>±, due to known astrophysical sources (at least for the 
electrons), and an unknown component, <j>±, in formulae: 
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the one needed to explain the 
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features in the spectra observed by PAMELA and 
FERMI-LAT. These experiments measure respectively 
the positron fraction and the total electron and positron 
fluxes as a function of the energy E of the detected , 
i.e.: 



P(E) = 



ME) 



<j>+{E) + <t>-{E) 



F(E) = <P+{E) + <f>-{E) . 

(2) 



The left-hand side of the equations above refer to the 
experimental measures. The contribution from the un- 
known source is then: 



cj>V(E) = P(E)F(E)-<P*(E) , 
4> U _{E) = F(E) (1-P(E))-^(E) 

In terms of their difference and sum: 
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= F(E) (2P(E)-1) + (^ B (E)-^(E)) , 



F{E)-(4> B {E) + ^{E)) 



(5) 



The latter equation implies F(E) > (j) B {E) + 4> B (E). 
We model the background spectrum using (j>±{E) = 
N B B ± (E), where N B is a normalization coefficient such 
that F(E)/(B-(E) + B+(E)) > N B and B ± (E) are pro- 
vided using specific astrophysical models. In this pa- 
per we adopt the popular Moskalenko and Strong model 
0[7], for which N B is less than 0.75 and B ± {E) are given 
by: 
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where 

GcV _1 cm _2 sec _1 sr _1 units. We checked that our results 
remain unchanged when replacing the parameterization 
above with the one adopted by the Fermi Collaboration 
(model zero) [ICTO]. 
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It is convenient to introduce the following parameter: 
<f>V(E) F{E)(l-P{E))-cj> B {E) 



ru(E) 



<t> u + {E) " P(E)F(E)-^(E) 



(10) 



This equation can be rewritten as 



R(E) 



F(E) l-(l + ru (E))P(E) 
B-{E) 



1 (v^l^ 



Nf 



(11) 



Although the sum rule R(E) seems to depend on the en- 
ergy it should, in fact, be a constant as is clear from the 
right hand side of the previous equation. This leads to a 
nontrivial constraint linking together in an explicit form 
the experimental results, the model of the backgrounds 
and the dependence on the energy of the unknown com- 
ponents. 

We now turn to the actual data and show in which way 
the sum rule (11) provides essential information on the 
unknown components of the CRs. Since we use simulta- 
neously the results of FERMI-LAT and PAMELA we are 
obliged to consider only the common energy range. Note 
that the CRs energy range below 10 — 20 GeVs, where the 
spectrum is affected by the Sun, is outside the common 
range. Within the relevant but limited range of energies 
we will consider here it is therefore sensible to assume rjj 
to be nearly constant. We find useful to plot the func- 
tion R(E) for different values of rjj in order to test the 
sum rule (111. This would imply that this function is in- 



dependent of the energy. The associated constant value 
would then be identified with the background CRs nor- 
malization factor Nb- We report the results in Fig. [T] 
The straight (red) line is the Nb = 0.75 value which is 
the largest one can assume for the background not to 
be larger than the FERMI-LAT results. We observe that 
there is a clear tendency for the combined data to predict 
a lower value of the constant Nb = 0.66,0.64,0.62,0.58 
for increasing value of the ratio rjj — 0,1/2,1,2. This 
is clear when looking, from top to bottom, at the dif- 
ferent panels of Fig. [T] It is interesting to note that we 
find a plateau, in the relevant energy range, up to rjj 
near the value of 2 when R(E) starts showing some devi- 
ation. The PAMELA and FERMI-LAT mean data have 
been used to determine the mean values of R(E) in the 
different panels of Fig. [T] As explained above we com- 
pare the data only in the energy range where the two 
experiments overlap. In the various panels we also show 
the uncertainties around the mean value which we have 
determined in the following way. We maximized (mini- 
mized) R(E) in (111 by using the one sigma deviations 



coming from both PAMELA and FERMI-LAT. This re- 
sulted in the shaded band around the mean value. Given 
the large uncertainties we cannot yet provide a more solid 
conclusion, however we can use the derived normalization 
for each different ratios of the unknown components to 
predict the positron fraction at energies higher than the 
ones provided so far by PAMELA. In order to be able to 
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FIG. 1: Ratio R(E) as a function of the energy E of electrons 
and positrons and for vales of ru = 0, 1/2, 1,2, from top to 
bottom. The shaded region accounts for the one sigma er- 
ror in PAMELA and FERMI-LAT. Secondaries are estimated 
according to the expressions in [5] [7]. 



make such a prediction we first rewrite ( 10 ) as follows: 

4> B (E) 



P(E) 



1 

l + r lr 



1 



F(E) 



(12) 



where we use for each rjj the estimated associated Nb- 
The different predictions for the positron fraction, as- 
suming that rjj remains constant over the entire energy 
range, up to 1000 GeV are shown in Fig. [2] The resulting 
picture disfavors both very small and very large values of 
ru which, in turn, means that one expects the electron 
fraction to neither be too small nor to indicate a large 
proportion of positrons. 
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FIG. 2: Model independent prediction for the positron frac- 
tion P(E) as a function of the energy E of electrons and 
positrons and for vales of ru — 0, 1/2, 1, 2, 4, from top to bot- 
tom. Secondaries are estimated according to the expressions 
in [HI [7] reported in the main text and the derived Nb values: 
N B = 0.66,0.64,0.62,0.58,0.50. 

Using the sum rule will allow to extract vital infor- 
mation from the data as they become more and more 
accurate. The special case rjj = 1 is an automatic pre- 
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FIG. 3: We display <f> u (E) = 4>±{E) corresponding to the 
central shaded region for the whole energy range of the 
FERMI-LAT experiment. The upper red points are the ac- 
tual FERMI-LAT data, and the background of electrons is the 
upper green stripe while the positrons is the lower magenta 
stripe. See the text for further explanation. 

diction of a great deal of models for dark matter which 
assume that charge symmetry holds both in the produc- 
tion and propagation of the unknown component of the 



CRs. 

It is, therefore, clear that our formalism is more uni- 
versal given that we have made no assumption in the 
derivation of the sum rules above. It is however, useful 
to adopt the model assumption of txj = 1 to derive the 
further constraint: 

me) . m^miim , (13) 

for which we now use N B = 0.62. We show in Fig. |] the 
sensitivity of the unknown flux (f> (E) = (j>±(E) on the 
experimental errors as well as on a 10% arbitrary varia- 
tion of the background spectrum. <\P (E) is the difference 
between the FERMI-LAT data (the red dotted points) 
and the positron (the magenta lower line) and electron 
(upper green line) backgrounds. The vertical spread as- 
sociated with the inner (blue) darkest region of <f> u (E) 
is due to FERMI-LAT error bars while the outer lighter 
shaded region also includes a 10% variation of the back- 
ground. This does not imply that the backgrounds are 
known to a 10% accuracy but merely shows how sensitive 
4> U (E) can be to such a variation. Although we investi- 
gated the case r u = 1 it is clear that our sum rules can 
be used to test any model of the unknown component. 

The sum rules introduced here are general and can 
be extended also to the proton and antiproton fraction. 
We have shown that current data still allow for approxi- 
mately equal contributions of the electrons and positrons 
from the unknown components of the associated CRs, but 
disfavor electron to positron fractions smaller than one 
half and larger than four. 

The current model independent analysis of the com- 
bined PAMELA and FERMI-LAT data shows that we 
will be able to deduce, thanks to the new constraints, 
vital information on the nature of the source and the 
propagation of the CRs. We have also demonstrated that 
the typical oversimplifying model assumption used so far 
constitutes a small portion of the allowed models still 
left unconstrained by the present data. Finally our sum 
rules can be easily used to constrain any specific model 
while we were able to predict, in a model independent 
way, the positron fraction at energies higher than the 
ones explored so far. 
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